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Abstract. The classical relationship between the Tutte polynomial of graph theory and 
the Potts model of statistical mechanics has resulted in valuable interactions between the 
disciplines. Unfortunately, it does not include the external magnetic fields that appear in 
most Potts model applications. Here we define the V-polynomial, which lifts the classical 
relationship between the Tutte polynomial and the zero field Potts model to encompass 
external magnetic fields. The V-polynomial generalizes Noble and Welsh's VF-polynomial, 
which extends the Tutte polynomial by incorporating vertex weights and adapting contrac- 
tion to accommodate them. We prove that the variable field Potts model partition function 
(with its many specializations) is an evaluation of the V-polynomial, and hence a polynomial 

i— i with deletion-contraction reduction and Fortuin-Kastcleyn type representation. This unifies 

an important segment of Potts model theory and brings previously successful combinatorial 

r\ machinery, including complexity results, to bear on a wider range of statistical mechanics 

models. 

This e-print is an extended version, including additional background information and 
more detailed proofs, of a paper of the same name that is to appear in Advances in Applied 
H Mathematics. 
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1. Introduction 
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The classical relationship between the Tutte polynomial of graph theory and the Potts 
model of statistical mechanics has resulted in valuable interactions between the fields, par- 
ticularly in the investigation of zeros and computational complexity. Unfortunately, the 
classical theory does not encompass the external magnetic field that appears in most Potts 
model applications. Our main result is a new graph polynomial, V, that lifts the well known 
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relation between the classical Tutte polynomial and the zero field Potts model to the full 
Potts model with external magnetic field, and provides a single polynomial for the study of 
previously disparate Potts models. It furthermore specializes to the previously known result 
relating the classical Tutte polynomial and the zero field Potts model. This provides a frame- 
work for extending results from Tutte polynomial investigations to the variable field Potts 
model, and vice versa. In particular, the relationship shows that the variable field partition 
function is a polynomial in a set of expressions, gives a deletion-contraction reduction that is 
conducive to induction arguments, and provides a Fortuin-Kasteleyn-type representation for 
the model. Existing results for the W- and [/-polynomials then make connections between 
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the variable field Potts model and knot theory, and also provide computational complexity 
results for variable field Potts models. 

The Potts model of statistical mechanics models how micro-scale nearest neighbour energy 
interactions in a complex system determine the macro-scale behaviours of the system. This 
model plays an important role in the theory of phase transitions and critical phenomena 
in physics, and has applications as widely varied as magnetism, tumor migration, foam 
behaviours, and social demographics. The underlying network of the system is typically 
modelled by a graph, with the vertices representing molecules, cells, bubbles, households, 
etc., with edges between 'neighbouring' units, i.e. those that may influence each other. The 
vertices may have different properties called spins assigned to them, representing anything 
from magnetic spin to socio-economic status of the household, and neighbouring vertices 
have an interaction energy. A state of a graph is a choice of spin at each vertex. The Potts 
model theory has its origins in the study of magnetism, and we use that terminology, for 
example spins and states, here. Two recent surveys of the Potts model and its connections 
with graph theory are [5] and [54|. 

The g-state Potts model partition function is the normalization factor for the Boltzmann 
probability distribution and is given by: Z(G) = Yl e - ^^)) where the sum is over all 
possible states a of a graph G using q spins. Here (3 is a function of the temperature, and 
the Hamiltonian h is a measure of the energy of the state. The Hamiltonian may have many 
different forms depending on the specific application. In the absence of an external magnetic 
field and with constant interaction energies between adjacent vertices, the Hamiltonian of a 

state a is simply h(a) = —J ^2 <5(ctj,ctj), where 5 is the Kronecker 5-function. In this 

{i,j}eE(G) 
case, the Potts model partition function is equivalent to the classical Tutte polynomial: 

Z (G; q, (3) = q KG) v HG)\-KG) T f G . 9±^ v + x 

where v = e J/3 — 1. (See jTTj for the nascent stages of this theory and later exposition in 
[U 13 13 El E3], EH E31 [M] ) • This relationship has resulted in a remarkable synergy between the 
areas of combinatorics and statistical mechanics, particularly for computational complexity 
and the study of the zeros of these polynomials. (See [5j [191 H21 El [55] for surveys of these 
results.) 

The equivalence of the Tutte polynomial and the Potts model partition function assumes 
the absence of an external field. However, many applications of the Potts model depend 
on additional terms in the Hamiltonian corresponding to the presence of such additional 
influences, for example the standard models of magnetism, the cellular Potts model of [18] . 
and also see [40J for examples in the life sciences. Furthermore, within many of these models, 
there is a need for edge dependent interaction energies and site dependent external fields. Sig- 
nificant work has been done in incorporating the edge dependent interaction energies, both 
from a combinatorial perspective with edge weighted generalizations of the Tutte polynomial 
(see [HI [151 [57]) an d from a statistical mechanics approach (see [37] for a survey). Only very 
recently though have external fields been investigated in the context of graph polynomials. 
In the special case of the Ising model, essentially the Potts model with q = 2, an interesting 
new polynomial has been found in [52] . This polynomial has a deletion-contraction relation 
and captures approximating functions, but does not specialize to the Tutte polynomial, nor 
does the deletion-contraction extend to non-constant magnetic fields. Also, some weighting 
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strategies have been developed to study non-zero magnetic fields in a graph colouring frame- 
work (see [lOj [HJ SI]), with an accompanying graph polynomial. These papers show that 
the Potts model partition function with an external field term does not have a traditional 
deletion-contraction reduction, and in fact give details for its deviation from this. 

Here however, by using a definition of contraction that incorporates vertex weights, we 
are able to assimilate a Hamiltonian of the following generic form into the theory of the 
Tutte-Potts connection: 

9 

(1) h(o) = - ^2 ■Jij$(°'ii cr j) - ^2 S ^2 M i, a 5{a,a i ) 1 

{i,j}£E(G) VidV(G) a=l 

where a magnetic field vector 

M i :={M i ,i,M <i2 ,...,M i> J 

is associated to each vertex V{. From this generic form we are able to specialize to various 
forms of the Hamiltonian with external fields that are common in the physics literature. 

Our approach is motivated by the W- and [/-polynomials of Noble and Welsh [37] . These 
graph polynomials were originally developed in the context of knot theory through an in- 
vestigation of the combinatorics behind Chmutov, Duzhin and Lando's work on Vassiliev 
invariants (see [T21 [T3"t IT^] ) . They generalize the classical Tutte polynomial by incorporating 
vertex weights, while retaining essential properties such a deletion-contraction reduction, 
state sum formulation, and universality (or recipe theorem). Like the classical Tutte poly- 
nomial, the jy-polynomial has a deletion-contraction reduction, but it begins with positive 
integer weights on the vertices, and these weights are summed if an edge joining them is con- 
tracted. We define a V-polynomial that extends the ly-polynomial by incorporating edge 
weights to encode variable interaction energies and allowing vertex weights in a semigroup, 
in particular the vector space C q . 

If we think of the weighted spins in the Hamiltonian as vertex weights, then the additional 
term ^2 i&v r G ) Yla=i ^i,a^{ a ^ a i) i s handled in the deletion-contraction relations by summing 
the vector valued vertex weights exactly as done by the U-, V-, and H^-polynomials. We 
use this insight about the vertex weights to prove that the Potts model partition func- 
tion with an external field is an evaluation of the V-polynomial. This gives the desired 
deletion-contraction reduction for the external field Potts model, facilitating induction ar- 
guments. More importantly, these partition functions may now be expressed as polynomials 
and Fortuin-Kasteleyn-type representations, which facilitates Taylor expansions and com- 
puter simulations used to estimate critical exponents and phase transitions. 

Since the classical Tutte polynomial is a specialization of the V-polynomial, this relation- 
ship contains the original relationship between the classical Tutte polynomial and zero field 
Potts model. This new relationship for the variable field Potts model extends the compu- 
tational and analytic tools available to statistical mechanics applications. In particular, we 
are able to immediately transfer computational complexity results for the W- and U- poly- 
nomials to the extended Potts model partition functions. Other results, such as those on 
zeros and phase transitions in the classical and multivariate settings, might also be adapted 
to this new context. Furthermore, provocative new connections between knot theory and 
statistical mechanics arise from this theory. 
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2. Background 

2.1. The g-state Potts model. We begin by recalling some basic information about the 
Potts model. Let G be a graph and consider a set of q elements, called spins. In the abstract, 
the spins may be numbers or colours, but typically they are values relevant to some specific 
application. For example, in studying uniaxial magnetic materials, q — 2, and the possible 
spins are +1 and — 1. A state of a graph G is an assignment of a single spin to each vertex 
of the graph. Combinatorially, a graph state may be thought of simply as a colouring (not 
necessarily proper) of the graph. Thus, if the vertex set of G is V(G) = {i>i, . . . ,v n }, then a 
state of G is an assignment a : V(G) — ¥ {1, . . . , q}, for q G Z + . We let S(G) denote the set 
of states of G. 

For greater readability we let <7j := <t(i>j), for er G 5(G). In addition, to simplify some of 
the summation formulas, at times we will use the indices % — 1, . . . , n of the vertices in place 
of the vertices, for example we may denote an edge e = {vi,Vj} by e = {i, j}, and a vertex 
V{ as i. This convention will always be clear from context and should cause no confusion. 

The interaction energy may be thought of simply as a weight on edge of the graph. In 
physics applications, the interaction energies are typically real numbers, with the model 
called ferromagnetic if all are positive, and antiferromagnetic if all are negative. However, in 
our context of graph polynomials, they may simply be taken to be independent commuting 
variables. We will denote the interaction energy on an edge e = {v j, Vj} by J e := Jij(— J Vi ,v), 
and let 

3W) : = - ^2 Jijbi.or^CTj) = - ^ JeKViiVj)- 

{i,j}&E(G) e£E(G) 

The Hamiltonian is a measure of the energy of a state. We begin with the simplest 
formulation, where the interaction energy is a constant J, and the Hamiltonian depends 
only on the nearest neighbour interactions (without any external field or other modifying 
forces) . 

Definition 2.1. The zero field Hamiltonian is 

(2) h(a) = -J J2 KWi), 

ijeE(G) 

where a is a state of a graph G, where o~i is the spin at vertex i, and where 5 is the Kronecker 
delta function. 

However, to encompass Hamiltonians with various other external fields and variable inter- 
action energies, we need a much more general form. 

Definition 2.2. Let G be a graph. Assign to each edge e an interaction energy of J e , and 
assign to each vertex Vi a magnetic weight vector 

M i :=(M itl ,M i , 2 ,...,M i , q ). 

Then the Hamiltonian of the Potts model with variable edge interaction energy and variable 
magnetic field is 

q 

(3) h(a) = - ^2 J i,jS( a ii a j) - ^2 ^2 M i,aS(a,ai). 

{i,j}£E VieV{G) a=l 
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Note that the Hamiltonian of Equation [2] may be recovered from the general Hamiltonian 
of Equation [3] by setting J e = J for all edges e, and taking all the magnetic field vectors 
to be zero. We will recover various other Hamiltonians of relevance to physics applications 
from the general Hamiltonian of Equation [3] in a similar way. 

Regardless of the choice of Hamiltonian, the Potts model partition function is the sum 
over all possible states of an exponential function of the Hamiltonian, as follows: 

Definition 2.3. Given a set of q spins and a Hamiltonian h, then the q-state Potts model 
partition function of a graph G is 

Z(G)= J2 ^^ 

ff€S{G) 

where /3 = 1/(kT), where T is the temperature of the system, and where k = 1.38 x 10~ 23 
joules/Kelvin is the Boltzmann constant. 

The Potts model partition function is the normalization factor for the Boltzmann prob- 
ability distribution. That is, the probability of a system in thermal equilibrium with its 
environment being in a particular state a at temperature T is: 

(4) Pr(a,P) = ex P (-ph(a))/Z(G). 

Important thermodynamic functions such as internal energy, specific heat, entropy, and 
free energy may all be expressed in terms of the partition function. An important goal of 
statistical mechanics is to determine phase transition temperatures, that is, critical tempera- 
tures around which a small change in temperature results in an abrupt, non-analytic change 
in various physical properties. Since the thermodynamic functions can be given in terms of 
the partition function, determining the analyticity of the partition function in the infinite 
volume limit (over a class of well-behaved graphs of increasing size), both theoretically and 
through computer simulations, is central to this theory. This has led to considerable study 
of the zeros of the partition function and its computational complexity. 

2.2. The U- and H^-polynomials. Since our methods build on Noble and Welsh's paper 
[3*7] , we recall some definitions and results from their work. 

As usual, E(G) and V(G) are the edge set and vertex set, respectively, of a graph G. 
If A C E(G), then n(A), r(A), and k(A) are, respectively, the nullity, rank, and number 
of components of the spanning subgraph of G with edges in A. A vertex weighted graph 
consists of a graph G, with vertex set V(G) = {vi,V2, ■ ■ ■ ,v n } equipped with a weight 
function u : V(G) — > Z + . The weight of the vertex Vi is the value u(vi). For simplicity of 
exposition we will often write u>i rather than U)(vi), and x for an indexed set of variables, 
e.g. x = {xk\k G S} for some set S. 

Definition 2.4. If G is a vertex weighted graph with weight function u, and e is an edge of 
G, then: 

(1) If e is any edge of G, then G — e is the graph obtained from G by deleting the edge 
e and leaving the weight function unchanged; 

(2) If e is any non-loop edge of G, then G/e is the graph obtained from G by contracting 
the edge e and changing the weight function as follows: if Vi and Vj are the vertices 
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Figure 1. Deletion and contraction of an edge in a vertex weighted graph. 

incident to e, and v is the vertex of G/e created by the contraction, then u(v ) = 
u(vi) + u{vj). Loops are not contracted. 

The deletion and contraction of an edge in a vertex weighted graph is illustrated in Fig- 
ure [TJ 

We now recall Noble and Welsh's VF-polynomial. 

Definition 2.5 (Noble and Welsh [37j). Let G be a vertex weighted graph and {y, X\, x 2 , ■ ■ ■} 
be a set of commuting indeterminates. The W -polynomial of G, 

W(G) =W(G;x,y) eZ[{ Xi } i&+ ,y] 

is defined recursively by 

(1) if e is not a loop then 

W{G) = W{G - e) + W{G/e), 



where deletion and contraction act as in Definition 2.4 
(2) if e is a loop, then 

W(G)=yW(G-e); 



(3) if E m consists of m isolated vertices of weights u>\, . . . , u m , then 

in 

W(E m ) = l[: 



in 
i=l 

The VF-polynomial satisfies a recipe theorem (or universality theorem) that is analogous 
to the recipe of the Tutte polynomial Q41J). 

Theorem 2.6 (Noble and Welsh [37]). Let f be a function on vertex weighted graphs defined 
recursively by the following conditions: 

(1) if e is not a loop then 

f(G) = af(G-e) + bf(G/e). 

(2) If e is a loop, then 

f(G) = yf{G-e). 

(3) If E m consists of m isolated vertices of weights u)\, . . . ,u m , 

rn 

X u 

i=\ 

Then, if a and b are non-zero, 

ax y 



f{E m ) = \\. 



/(G) = a mG)nv(G)\ b mG)\ w ^ G . ^ 



(i 



For reference later, we record the following properties of the jy-polynomial. 

Theorem 2.7 (Noble and Welsh [37J). Let G be a vertex weighted graph. Then 

(1) W(G) can be represented as a sum over spanning subgraphs: 

W(G)= J2 x Cl x C2 ---x Ck{A) (y-ir^ A \ 

ACE(G) 

where q is the sum of the weights of all of the vertices of the l-th component of the 
spanning subgraph (V(G),A). 

(2) If Xi = 9, for each i, then the resulting polynomial W(xi = 8,y) is given by the Tutte 
polynomial: 

W(x i = e,y)=9 k ^T(G;l + 9,y). 

The U -polynomial, implicitly studied by Chmutov, Duzhin and Lando in [13] , was formally 
defined by Noble and Welsh in [37] as an adaptation of the IV-polynomial to unweighted 
graphs. The [/-polynomial of an unweighted graph is is defined simply by assigning the 
weight 1 to each vertex and then taking the jy-polynomial of the resulting weighted graph. 

3. The V-polynomial for edge weighted graphs 

In this section we define the V-polynomial, which extends the TV-polynomial to graphs 
that have both edge weights and more general vertex weights. An edge weighted graph G 
is a graph equipped with a mapping 7 from its edge set E(G), to a set 7 := {7 e } e e£(G)- 
Throughout this paper we will use the standard convention 7 : e 1— > j e , for e G E{G). We 
incorporate the edges weights using what has become a 'standard' technique (see for example 
[El |33l S3 EZ]) for extending Tutte-like polynomials to weighted graphs: we associate an 
indeterminate to each edge that then becomes a factor in the deletion-contraction reduction. 

We also generalize the vertex weights. In the original definition of the IV-polynomial (|_37J) 
the vertex weights are taken to be in Z + . Here however, we only require that the weights be 
in a torsion-free commutative semigroup (we have not explored the possibility of torsion, but 
this would be an interesting investigation). All that is actually necessary is that the vertex 
weights be additive and be a subset of the indexing set S of the variables x = {xk}kes of 
the polynomial. In fact, the weights need not range over the whole indexing set, provided 
that every possible sum of the vertex weights in each component of G appears as an index in 
x. For the physics applications later in the paper, we will usually consider vertex weights in 
a g-dimensional complex vector space. In the definition of V, as with W, when a non-loop 
edge is contracted, the vertex weights on its ends are added. 

The letter V is logical in the progression U-, V-, W-, and we use boldface simply to 
distinguish the polynomial V(G) from the vertex set V(G). 

Definition 3.1. Let S be a torsion-free commutative semigroup, let G be a graph equipped 
with vertex weights w := {oJi} C S and edge weights 7 := {7 e }, and let x = {x k } keS be a 
set of commuting variables. Then the V -polynomial of the vertex and edge weighted graph 
G, 

V(G) = V(G,w;a5,7) e Z[{7e}ee£(GQ,{zfeW], 

is defined recursively by: 



(1) if e is a non-loop edge of G, 

V(G) = V(G-e)+<y e V{G/e), 

where contraction acts as described in Definition |2.4| 

(2) if e is a loop, then 

V(G) = ( le + l)V(G-e); 

(3) if E m consists of m isolated vertices of weights u>\, . . . , u m , then 



V(E m ) = J] 



in 
.i u 

i=l 



We now prove several basic results about the V-polynomial. The proofs of most of these 
results are adaptations of Noble and Welsh's proofs of the corresponding results about the 
iy-polynomial from [37] . 

Proposition 3.2. The polynomial V is well defined in the the sense that the polynomial is 
independent of the order in which the deletion-contraction relation is applied to the edges. 

Proof. We proceed by induction on the number of edges of G. The result clearly holds for all 
graphs with fewer than two edges. Assume that G has at least two edges, and that e and / are 
distinct non-loop edges of G. First we calculate V(G) by applying the deletion-contraction 
relation to the edge e first, and / second: 

(5) V{G) = V(G-e)+ le V(G/e) 

= V([G-e]-f)+ V V([G - e]/f) + le V{[G/e] - f) + lelf V{[G/e}/ f). 

Applying the deletion-contraction relation to the edge / first, and e second gives 

(6) v(G) = V(G-f)+ lf V(G/f) 

= V([G-f]-e)+ le V{[G - f]/e) + V V([G/f] - e) + lfle V{[G/ f]/e) 

Then, since [G - e] - / ^ [G - f] - e, [G - e]/f = [G/f] - e, [G/e] - f = [G - f]/e and 
[G/e}/f = [G/f]/e, the right hand sides of ^ and ^ are equal. A similar argument shows 
the independence of order when e or / is a loop. The result then follows by induction. 

□ 

Theorem 3.3. V(G) can be represented as a sum over spanning subgraphs: 

' \.G) = / j X Cl X C2 ■ ■ ■ Xc k{A) Y J_ 7e> 
ACE(G) eeA 

where c\ is the sum of the weights of all of the vertices in the l-th connected component of 
the spanning subgraph (V(G),A). 

Proof. Let 

S\G) '■= 2_^t X ci X c 2 " ' X c k(A) J_J_7e- 
ACE(G) eeA 

We will show by induction on the number of edges of G that S(G) = V(G). If G has no 
edges then the result is clearly true. So suppose that e is a non-loop edge of G. Then 
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(7) S(G) = Yl x <* x >* ■ ' ' X ^(A) II 7e 

ACE(G) egA 

\ / 

/ j x c 1 x c 2 ' ' ' x c k ( A) J_ J_ 7e + / j x a x c 2 •■• %c k ( A ) J_J_ 7e 



ACB(G) eeA 

e£A 



/ 



i ^ce(G) eeA / 

\ eeA / 



There is a natural bijection between the spanning subgraphs of G that contain the edge e 
and the spanning subgraphs of G/e (given by A K A - { e })- Using this correspondence 
and the obvious correspondence between the spanning subgraphs of G — e and the spanning 
subgraphs of G that do not contain e, we can write the above expression as 

/ , X ci x c 2 ' ' ' X c k ( A ) J_ J_ 7e I + 7e I / j x c 1 x c 2 ' ' ' X c k ( A) J_ 1 7e 

K ACE(G-e) eeA J \ACE(G/e) eeA 

This is equal to S(G — e)+j e S(G/e). A similar argument shows that S(G) = ( , y e + l)S(G — e) 
when e is a loop. It then follows by induction that S(G) = V(G). □ 

We now give a recipe theorem for V. 

Theorem 3.4. Let f be a function on vertex and edge weighted graphs defined recursively 
by the following conditions (where each a e ^ 0): 

(1) for any non-loop edge e 

f(G) = a e f(G-e) + PJ(G/e), 

with a e 7^ for all e; 

(2) for any loop e 

/(G) = (a e + &)/(G-e); 

(3) If E m consists of m isolated vertices of weights ui±, . . . ,oo m , 

m 

X„ 
i=l 

Then 



f(E m )=Y[: 



f( G )=[ 11 a e\V(G,u;x,{l3 e /a e } e€E{G) ) 

\e&E(G) J 

Proof. Let 

/:= l /, 

(UeeE(G) a e) 

and let 7 e := (3 e /a e . Then 

f(G) = f(G-e) +1 J(G/e), 
when e is a non-loop edge of G, and 

7(G) = (l + 7e )7(G-e), 

when e is a loop. So / is equal to the V-polynomial and the result follows. □ 
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The multivariate Tutte polynomial (see [37J and also [3H]), which we will denote by Zt to 
distinguish it from the partition function Z, is an extension of the Tutte polynomial to edge 
weighted graphs, defined by 

Z T (G;9,j): = Y, O k{A) \[le- 

ACE{G) eeA 

In the following theorem we prove that the V-polynomial generalizes the IV-polynomial 
and the multivariate Tutte polynomial. We first show that the VT-polynomial can be re- 
covered from the V-polynomial by setting all the 7, variables in the V-polynomial equal to 
y — 1. Then, given that the V-polynomial is an edge weighted version of the P^-polynomial, 
and that the VT-polynomial extends the classical Tutte polynomial, it is hardly surprising 
that the V-polynomial also extends the multivariate Tutte polynomial. 

Theorem 3.5. The V -polynomial generalizes both the W -polynomial and the multivariate 
Tutte polynomial. 

(1) If co : V(G) -> Z+ and j e = (y - 1) for each e G E{G), then 

V(G,oj;x, le = (y- 1)) = (y-lf^W (g,u; ^,v) ■ 

(2) If X{ = 6, for each i G Z + ; then, independent of the weights 00, 

V(G,u; Xi = 6,~f) = Z T (G;9,~f). 



Proof. For the first item, we apply the recipe theorem for V (Theorem 3.4) to 



f{G) = {y-l)^c)\ W ( G .^ y 



x 



By Definition 2.5 



(1) f(G) = f(G -e) + {y- l)f(G/e) for e a non-loop edge; 

(2) f(G) = yf{G - e) if e is a loop; 

(3) /(£„) = niLi^. 



Thus, / satisfies the criteria of Theorem |3.4| with a e — 1 and (3 e = y — 1, for all e, and the 
result follows. 

The second item follows easily from the state sum for the V-polynomial: 

V(G;x, t = 9 n )= J2 fll e ](n^) = E O k{A) \{le = Z T {G-e n ). 

ACE(G) \i=l J \e&A J A<ZE{G) eeA 



The second item also follows from the first item and the second item in Theorem | 

□ 



Although Theorem 3.5 tells us that the V-polynomial assimilates the multivariate Tutte 
polynomial Zt, the relation between the V-polynomial and the fully parameterized coloured 
Tutte polynomial (see [HI HHJ [57] ) remains to be explored. 
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4. The V-polynomial and the Potts model in an external field with edge 

dependent interaction energies. 

We now come to our reason for creating the V-polynomial. Here we show that the Potts 



model partition function with the Hamiltonian h(a) from Definition 2.2 is an evaluation of 
the V-polynomial. This provides a recursive deletion-contraction definition for the full Potts 
partition function with variable edge interaction energy and variable magnetic field, and also 
shows that it is a polynomial. 

Theorem 4.1. Let G be a graph equipped with a magnetic field vector M^ = (M^i, . . . , M itQ ) G 
C 9 at each vertex Vi, and suppose 

g 

h{a) = - ^ Ji,j S ( a h^j) ~ ^2 ^2 M i,a$( a i a i)- 
{i,j}€E Vi£V(G) a=l 

Then 

(1) if e = {vi,Vj} is a non-loop edge ofG, 

Z[G) = Z[G - e) + (e^ - l)Z{G/e); 

(2) if e = {vi,Vi} is a loop 

Z(G) = e^Z{G-e); 

(3) Z(E n ) = Yl^i-X-Mi, where E n consists of n isolated vertices of vector valued weights 

Mi, M 2 i • • • , M n , and, for any weight M it 

<? 

pM i>a 



a=l 

Proof. To prove Item HI of the theorem, let e = {v a ,Vb} be a non-loop edge of G. Then 



Z(G)= Y^ e~ ma) = Y^ e -Ptt°) TT qE„ £»*,«*(«.<*) 

ctS5(G) cts5(G) WeV(G) 



E ***> n 

ffe-S(G) WeV(G) 



^Mi, . 



The last equality holds since Mi >a S(a, a*) = M^ ai if a = a iy and is zero otherwise. 

We can separate the above sum by collecting together the states in which a a ^ cr&, and in 
which a a = a& to get 



J2 e-^' (<T) 


n etsm ^ 


+ £ e-*> 


n * 


treS(G) 


\meV(G) , 


1 <t£S(G) 
0'a=(Tf > 


\viev(G) 



PM % , 
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If a is any state of G, then there is a unique state a' of G — e where each vertex has the 
same spin as it has in a. Using this correspondence, the expression above then becomes 



J2 e-^ (CT) 


n em ^ 


+ e^ J2 e '^ {a) 


17 e^*« 


a-eS(G-e) 


Wv(G-e) J 


' aeS(G-e) 


Wv(G-e) 


0"a7^0"(, 


\ / 


O a =0~b 


\ 



The left hand sum in the above expression is nearly Z(G — e), but we are missing the states 
of G — e where a a = at,. So we simply add and subtract these states, to obtain 



Z(G) 



£ 



-PKa) 



n 



JM it , 



+ e 



/3J e 



D £ 



-tfito 



n 



,/3M,, 



treS(G-e) 



.Vi£V(G-e) 



o-eS(G-e) 
CTa=0'b 



,^eV(G-e) 



The left-hand side is precisely Z(G — e). 

We now turn our attention to expressing the right-hand sum as Z(G/e). If a is any state 
of G — e with a a = cTfe, then there is a unique state a" of G/e, where the vertex resulting 
from identifying a and b has the common spin a a = at,, and each remaining vertex of G/e 
has the same spin as it does in a. We immediately have j(a) = j(u"). Now for a state a 
with a a — at,, we have 



/ 



n 



JM ii0 



^SV(G-e) 



n 



JM it , 



. »,£V(G-e) 

\ j^a or b 



\ 



J 



( 



J3(M a ,<r a +M bia ) 



n 



^M<, 



\ i^a or fe 



\ 



/ 



2 /3(M a , CTa +M 6 , CTa ) 



where the last equality holds because a a = &b- Notice that M atUa + M^ aa is the <r a -th term 
of the vector M a + Mb- If v c is the vertex of G/e resulting from identifying v a and Vf,, then 
v c has weight M" c = M a + Af &. Also, in the state a' of G/e induced by the state a of G 
as described above, the spin of v c is equal to the spins of v a and v^. Using the fact that 
c" c = c"fe = 0"a; so that M aj(7a + M bjCra = M C)(Tc , we can write the above product as 



n 



^Afi, 



Z(G/e). 



Vi£V(G/e) 

Equation [8] then gives 

Z{G) = Z(G - e) + (e^ - l)Z(G/e) 
as required. 

Now, to prove Item|2]of the theorem, suppose that e = {v a ,v a } is a loop. Then 



z( G ) = y: ^ 



{i,j}€E(G) Ji,3 $(< T i> CT i) 



<re5(G) 



n 

,viev(G) 



JMi. 



Since e is a loop, both its ends are assigned the same spin in each state, so we can write the 
partition function as 



JJ, 



aeS(G-e) 



/E{ij}gB(G-e) J i,3 S i,3 



n 



JM t . 



.Vi&V(G-e) 



12 



where we have identified the states of G and G — e in the way described in the proof of 
Item u\ above. Thus, Z(G) = eP Ja < a Z(G — e) as required. 



Finally, if G consists of a single isolated vertex, it is easily checked that Z(G) = Yla=i e^ Ml ' a . 
Item (J3l) then follows by the multiplicativity of the partition function. □ 

The following theorem is the main result of this paper. It states that the Potts model 
partition function with a variable external magnetic field and variable edge interaction is an 
evaluation of the V-polynomial. 

Theorem 4.2. Let G be a graph equipped with a magnetic field vector Mi = (M^i, . . . , M itQ ) G 
C 9 at each vertex Vi, and suppose 

h{cr) = - ^ JijtifaiVj)- Y ^2 M i,aS(a,ai). 

{i,j}£E VieV(G) a=l 

Then 

Z(G) = V (G,w; {X M } MeCq , {e^- - % J}e£(G) ) , 

where the vertex weights are given by u(vi) = Mi and, for any M = (Mi, . . . , M q ) e C q , 

q 



x M - Yl 



e /3M a 



a=l 



Proof. The equality of the Z(G) and the V-polynomial follows immediately from Theo- 
rem 4.1 and the recipe theorem for the V-polynomial (Theorem 3.4). □ 



As a corollary to this theorem we obtain a Fortuin-Kasteleyn-type representation for the 
Potts model with variable external magnetic field and variable edge interaction. 

Corollary 4.3. Let G be a graph equipped with a magnetic field vector Mi = (M^i, . . . , M^ q ) e 
C 9 at each vertex v i} 

Q 

h(cr) = - Y Jij5{ou°j)- JZ ^2 M i,a$( a i a i)- 

{i,j}£E Vi£V(G) a=l 



Then 



Z(G) = £ X Mci ■ ■ ■ X MCk(A) ( l[(e^ - 1) ] 

ACE(G) \e&A J 



ACE{G) 

where Mc l is the sum of the weights, M i; of all of the vertices Vi in the l-th connected 
component of the spanning subgraph (V(G), A), and Xm = J2 9 a =i e /3Ma . 



Proof. The result follows immediately from Theorem |4.2| and Theorem |3.3| D 



Since the V-polynomial of a graph G is a polynomial in Z[{xjt}, {7 e }], we can use the 
relation between the V-polynomial and the Potts model to express the g-spin Potts model 
partition function with variable edge interaction energy and variable magnetic field as a 
polynomial. 
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Corollary 4.4. For a connected graph G equipped with a magnetic field vector Mi = 
(Mi : i, . . . , Mi >q ) G C q at each vertex Vi, the q-spin Potts partition function with variable 
edge interaction energy and variable magnetic field is a polynomial in the variables 

{v,X M \MeM}, 

where v = {e 0Je - l} eeE (G), and X M = Yl=i ePMa f or 
M:= {El=i G)l M*i|ei = Oorl}. 

This result may be extended to graphs with more than one component by using the 
multiplicativity of V. 

5. Some specializations 

As noted in the introduction, the V-polynomial provides a unifying framework for Potts 
models with various Hamiltonians. Here we highlight some special forms of the magnetic 
weight vectors and resulting Hamiltonian that are frequently studied in the physics litera- 
ture. We show how the associated partition functions are evaluations of the V-polynomial 
and its specializations. This implies that these important models are also polynomials with 
deletion-contraction reductions and Fortuin-Kastelan-type representations. We also show 



how Theorem 4^ can be applied to the various Hamiltonians to give a hierarchy of rela- 
tions among the resulting partition functions and the W-, U-, multivariate Tutte, or Tutte 
polynomials. These graph polynomials are all specializations of the V-polynomial. 

5.1. A hierarchy of specializations. The following result says that if we restrict to the 
Potts model with constant interaction energies and where all the magnetic field vectors are 
positive integer multiples of a given magnetic field vector, then the partition function is 
described by the VF-polynomial. 

Theorem 5.1. Let B G C q , and let G be a graph equipped with a magnetic field vector 
Mi := kiB = ki(Bi, . . . , B q ) at each vertex Vi, where ki G Z + . If the Hamiltonian is 

h{a) = -J 22 ^( a ^ a i)~ ^2 S 22 k iB a 5(a,a i ), 

{i,j}<=E(G) VieV(G) a=l 



then 
where 



Z(G) = (e^ J - iy v ^W (G>; {x k } k&+1 e^ 7 ) , 






!3kB a 



and B a is the ct-entry of B. 



Proof. Setting each Jj ,• = J in Theorem 4.2 gives 



Z{G) = V (G>; {A M } MeZ+ ,, 7e = e? J - l) 



where the Xm are defined in Theorem |4. 2 Each index M , of Xm , is then of the form kB 



for some k G Z + . Therefore Xm = X k B = X]'=i e^ kBa , where B a is the a-entry of B. 
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Also since every index is of the form kB, the indexing set is isomorphic to Z + . Thu s, we 
can choose our indexing set to be Z + , letting X k = X]«=i e /3fcBa . Applying Theorem 3.5 then 
gives 

Z (G) = (e? J - 1)\ V ^\W (G>; {x k } k ^ + , e? J ) , 

where x k = (^Li ePkBa ) I {^ J ~ 1 )- D 

If, in addition to constant edge interaction energy, all of the magnetic field vectors are 
fixed, then the Potts model partition function can be recovered from the ^/-polynomial. 

Theorem 5.2. Let G be a graph equipped with a fixed magnetic field vector B = (B\, . . . , B a ) 
at each vertex, and Hamiltonian 

<? 
h(a) = -J ^2 $(°~i,°~j)- Yl ^B a 5(a } o-i) } 

{i,j}GE(G) ViEV(G) a=l 

then 

Z (G) = (e" J - l) |y(G % (G>; {x k } k&+ , e? J - l) , 
where 



hi 






q 

j3kB a 



and B a is the a-entry of B. 



Proof. The result follows easily from Theorem 5^ upon noting that the vector IB is the 
vertex weight of each vertex of G, and recalling that initial weights of 1 on all vertices is the 
defining property of U . 

□ 

The well known classical results which relate the Potts model partition function and the 
Tutte or multivariate Tutte polynomial in the cases when M = may also be recovered from 
Theorem 14.11 and Theorem 13.51 as follows. 

Theorem 5.3. Let G be a graph. 

(1) Ifh{a) = -Y,{U\*E J iJ 5 (?i>< r 3)> th 



len 



Z{G) = Z T {G;q,{eP J u - l} {itj}eE(G) ); 
(2) and if h(a) = -JY.{i,j}eE S ( a ^ (J j)' then 

Z(G) = q k ( G ) v \ v (G)\-k(G) T ( G . ( q + v y v ^ v + 1 ^ 

where v = e^ J — 1. 



From the preceding theorems, we have the following hierarchy of relations among the V-, 
W-, U-, and Tutte polynomials, and the partition functions with Hamiltonians of various 
degrees of generality. This hierarchy holds in the case of the indiscriminately preferred 
spins given in Theorems 5^ through |5. 3 and adapts to the single preferred spin models of 



Theorems |5.4| through |5\6| presumably with other models following a similar pattern as well. 
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Z with variable interaction c - 
energies and zero field 




3 Z with variable interaction 
energies and variable field 



3 Z with constant interaction 
energies and integer scaled field 



--' Z with constant interaction 
energies and constant field 



3 Z with constant interaction 
energies and zero field 



5.2. Some more examples. There many variations of the Hamiltonian, having far too 
many forms to address here. Thus, we select the following frequently studied forms of the 
Hamiltonian simply to illustrate how the theory presented here may be adapted to other 
settings. 

One example of a Hamiltonian that is of particular interest in the physics literature mod- 
els a system with an external field in which one particular spin (the first, without loss of 
generality) is preferred (see for example the surveys [6], [55]). Here again, variable interaction 
energies and magnetic field values are allowed. In this case we have the following relation 
with the V-polynomial. 

Theorem 5.4. Suppose a complex value Z{ is associated to each vertex v-i of a graph G, and 
the Hamiltonian is given by 



K°) = ~ 7, Ji,j5(o-i, o-j) - 2_j z ^( 1 ' 

{i,j}£E Vi£V(G) 



o, 



Then 



Z(G) = V (G,w; {X z } zeC , {e 0Je - 1W (G) ) , 
where the vertex weights are given by u(vi) = Zi and 

X z = e? z + q- 1. 

Proof. We associate a magnetic field vector, Mi = (zi, 0, . . . , 0), with each vertex Vi of G. If 
we set 

h(a) := - ^ Jij&faiVj)- ^2 ^2 M i,aS(a,a,i), 

{i,j}&E(G) v^V(G) «=1 
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then 

Z(G)= J2 e~^ = Y. e " MW - 

aeS{G) aeS(G) 



Applying Theorem 4.2 with the Hamiltonian h yields that 



Z{G) = V (G>; {X M } MeCq , {e^ - l} e6S(G) ) , 
where the vertex weights are given by u)(vi) = M{ and, for any M £ C 9 , 



A M - ^ 



9 



C 



where M a is the a-entry of M. However, the only indices M £ C q that actually appear are 
of the form {z, 0, ... , 0) for some 2 £ C, and if TVf has this form, then Xm = Q l3z + Q — 1. 
Thus, we can choose our indexing set to be C instead of C 9 , so that 

Z(G) = V (Ceo; {X z } zeC , {e? J ' - l} eeS(G) ) , 

where X z = e^ z + q — 1, as claimed. □ 



An analog of Theorem 5.2 and associated hierarchy may easily be derived: if all of the z^'s 
are equal and Jjj = J, then Z(G) can be recovered from the [/-polynomial; and if all of the 
Zi's are of the form kiZ, for some fixed z £ C with the fcj's in Z + , and Jjj = J, then Z(G) 
can be recovered from the H^-polynomial. 

When the Hamiltonian is that of Theorem 5.4 a Fortuin-Kasteleyn-type representation 
for the Potts model is well known. This Fortuin-Kasteleyn-type representation, given in 



Corollary 5.5, can be immediately recovered from Corollary 4.3 



Corollary 5.5. Suppose a complex value Zi is associated to each vertex Vi of a graph G, and 
the Hamiltonian is given by 

h{cr) = - 22 JijtiivhVj)- X^ z iKh°i)- 

{i,j}eE viev(G) 

Then 

A<ZE{G) VeeA 

where Zc x is the sum of the weights, Zi, of all of the vertices Vi in the l-th connected component 
of the spanning subgraph (V(G), A), and X z = q — 1 + e^ 2 . 



In [45], Sokal studied the Potts model partition function in the case where the magnetic 
field vectors are of the form Mi = (M^i, Mj )2 , . . . , M^ r , 0, . . . , 0), where < r < q is fixed. 
With these magnetic field vectors, Theorem 4^ immediately gives the following result. 

Theorem 5.6. Let G be a graph where each vertex Vi is equipped with a magnetic field vector 
Mi = (M i}1 , M i>2 , • • • , M i>r , 0, . . . , 0), and suppose 

h(cr) = - Y Ji,j5(o-i,o-j) - Y ^2M iia 5(a,(Ti). 

{i,j}eE(G) v z eV(G) a=l 
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Then 

Z(G) = V (G,u; {X M } MeCq , {e^ - l} {a , b} mG)) , 
where the vertex weights are given by uj(vi) = Mi and, for any M e C q , 



X M = q-r + ^ 



e /3M a 



o=l 



where M a is the a-entry of M. 



Again, an analog of Theorem 5.2 and associated hierarchy for this external field may 
easily be derived: if each Mi = kiB, for ki G Z + , then Z(G) can be recovered from the 
VF-polynomial; and if all of the M{ are equal then Z(G) can be recovered from the U- 
polynomial. 

Also in [15] , Sokal found a Fortuin-Kasteleyn-type representation for the partition function 
used in the above theorem. This Fortuin-Kasteleyn-type representation can be immediately 
recovered from Corollary |4.3| 



Corollary 5.7. Let G be a graph where each vertex Vi is equipped with a magnetic field 
vector Mi = (M^i, M^, • • • , Mj jr , 0, . . . , 0). In addition let 

_ q 

h(o) = - 22 J i^{o-i,0-j) - 22 z2 M i,a^( a j a i)- 

{i,j}€E(G) ViEV(G) a=l 

Then 

z{ G )= y: x M Cl ---x MCk{A (n(e^-i) 

ACE(G) \eeA 

where Mc t is the sum of the weights, Mi, of all of the vertices Vi in the l-th connected 
component of the spanning subgraph (V(G), A), and Xm = q — r + Y2a=i e^ Ma . 

Our final examples are variations of the Ising model (essentially the q = 2 Potts model) 
used to study glassy behaviours (see [211 EEj)- The first is the Ising spin glass model, 
which has edge dependent random bond strengths J e , but no external field. The second is 
the Random Field Ising Model (RFIM), used to study disordered states. It has a random 
magnetic field in that the z^'s are randomly chosen local magnetic fields that each affect 
only a single site. To avoid redundancy of proof, we merge the two models in the following 
theorem (this generalization is also sometimes called the RFIM), and then recover each of 
them separately as corollaries. We choose these examples in part because complexity results 
for Ising models (or q = 2 Potts models) differ significantly from those for q > 2 Potts models, 
both in the classical (no external field) case and the examples here. Furthermore, there have 
been two very recent studies of the Ising model from a graph theoretical perspective. (3] treats 
the Ising model with constant interaction energies and constant (but non-zero) magnetic field 
as a graph invariant and explores graph theoretical properties encoded by it. [52] on the 
other hand, creates a new graph polynomial, which has a deletion-contraction reduction for 
non-loop edges, and which is, up to a prefactor and change of variables, equivalent to the 
RFIM. We discuss this polynomial further in the conclusion. 

The Ising model takes spin values in { — 1,1}. We will let r denote a state for the Ising 
model, which is a map r : V(G) — > { — 1,+1}. As usual we set Tj := r(vi). Also we will let 
T(G) be the set of states for the Ising model. 
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Theorem 5.8. Let G be a graph with a vertex weight u(vi) = Zi G C associated to each 
vertex Vi, and suppose the Hamiltonian and partition function are given by 

(9) h(r) = - J2 JijTiTj- J2 Z ^ and Z ^= J2 e ~ mT) - 

{i,j}€E(G) Vi£V(G) reT(G) 

Then 

Z{G) = e-K^^+sE^c^V (G,u; {x z } zeC , {e 2 ^ - l} eeB(G) ) , 

where for any z G C, x z = e 2z + e 4z . 

Furthermore, Z{G) is a polynomial in the variables {v,x z \ z G M.} , where v = {e l3Je — 



l}ee£(G); and M = |^!=i G)l £ i z i 



£i — or 1 >. 



Z(G) also has Fortuin-Kasteleyn-type representation 

ACE(G) \eeA 

where zq is the sum of the weights, z^, of all of the vertices Vi in the l-th connected component 
of the spanning subgraph (V(G),A). 

Proof. There is a natural bijection between T(G) and S(G) = {a : V{G) — > {1,2}}, which, 
given r G T(G), is determined by setting o~i = 1 if r, = —1, and Oi = 2 if r» = +1. This 
determines a state a G 5(G). 

Under this bijection between T(G) and S(G), observe that Ti = 1o~i — 3 and t^Tj = 
2<5(<7j, o~j) — 1, and so we may write the Hamiltonian in Equation fcty as 

h(cr) = -2 ^J Ji,j5(ai, o-j) - 2 ^ ^a» + ^ J»j + 3 ^ z* 
{i,j}sB(G0 KeV(G) {i,j}&E ieV(G) 

With this Hamiltonian, Z(G) = E re r(G) e_/3MT) = E ct& s ( g) e~^ ((T) . 
We now let 

{i,j}£E i&V{G) 

Then 

Z(G) = e _/3 (^ eei5 Je+3 E lG v(G) 2: y^ e -/ ^ (CT) . 

<re.S(G) 
By applying Theorem |4.2 we then have 

Z(G) = e-^«- ^^ £ v (G )^)y ( GjW . {X M } MeC2 , {e 2 ^ - l} eeB(G) ) , 

where the weight function u; is now Wj = A4"j := 2^(1,2) G C 2 , and, if M = z(l,2) then 

Finally, observing that since each M = z(l,2), for some z G C, we can take C as the 
indexing set instead of C 2 . This gives, 

Z(G) = e<^ J « + ^v«» Zi )v (G>; {x z } z£C , {e 2 ^ - l} ee£(G) ) , 

where x z = e 2z + e Az . This gives the relation between the V-polynomial and Z(G). 

The proofs that Z(G) is a polynomial in the given expressions and has a Fortuin-Kasteleyn- 
type representation are also very similar to those given previously, and are therefore omitted. 

□ 
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We first recover the well known result for the spin glass model. 
Corollary 5.9. The Ising spin glass model given by 

h(a) = - Y, J i,3 T i T j> and z (°) = J2 e_/3/l(r) ' 

{i,j}£E(G) reT(G) 

has partition function 

Z{G) = e-^- J *Z T (G; q = 2, {e 2 ^ - l} eeE{G) ). 
More importantly, we have the following results for the RFIM. 
Corollary 5.10. The RFIM given by 

h(r) = -J J2 r * r i ~ J2 ZiTi > and Z ^ = Yl e ~ mT) > 

{i,j}eE iaV{G) reT{G) 

has partition function 

Z(G) = e -«W)l+ 3 ^c«)^y (G,u- {x z } zeC , e 2/3J - l) , 

where the vertex weights u are given by u(vi) = Zi, and, for any z G C, x z = e 2z + e 4z . 
Furthermore, Z{G) is a polynomial in the variables je' 3,7 — l,x z | z G Aij , where M. = 

{Ell[ G)l ^\e t = 0orl}. 

Z(G) also has Fortuin-Kasteleyn-type representation 

Z(G)= J2 ^•••^o fc(A) (e /3J -l) |A| , 

ACE(G) 

where Zq is the sum of the weights, Zi, of all of the vertices Vi in the l-th connected component 
of the spanning subgraph (V(G),A). 

6. Leveraging Prior Results 

6.1. Computational complexity. Realizing the Potts model partition function as an eval- 
uation of the V-polynomial now means that the computation complexity results for the W- 
polynomial apply directly to partition functions with an external field. We collect some of 
these here, drawing on the work of Noble and Welsh [37], and Noble [3"8] . 

The computation complexity of the Tutte polynomial, and hence zero field Potts model 
partition function, has been extensively studied (see, in particular, Jerrum [23], Jaeger, 
Vertigan, and Welsh [22] and also Vertigan [5T]). The conclusion of these investigations 
is that computing the Tutte polynomial is #P-Complete for general graphs, except when 
q = 1, or for a small set of special points, or for planar graphs when q = 2. Furthermore, 
approximation is provably difficult as well: see [TJ El [M] for overviews and [19j [20], El] for 
recent results in this area. 



Thus, it is not surprising that the computational complexity consequences of Theorem 4.1 
and Section [5] for the variable field Potts model are somewhat bleak. Noble and Welsh [3Z 
have shown that computing any coefficient of the jy-polynomial is (jP-hard even for trees, and 
specific coefficients are (jP-hard for complete graphs. Thus the complexity of the variable field 
Potts model is at least as problematic (presumably more so if variable interaction energies 
are also used). Additionally, Noble and Welsh [37] have shown that computing evaluations 
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of the VT-polynomial, and hence the Potts model partition function with external field, are 
JjP-hard not only for trees, but even just for stars. 

The prognosis in the case of a constant, but non-zero, magnetic field, is somewhat better. 
As noted above, the ^-polynomial corresponds to a constant magnetic field vector. Noble 
[38] has shown that if G is a graph with tree-width at most K, then the [/-polynomial, and 
hence the partition function of G, may be evaluated at any point in roughly 0(a>Kn 2K+3 ) 
arithmetic operations. 

Complexity results for the Ising model, essentially the q = 2 Potts model, differ signif- 
icantly from the general Potts model. In particular, the partition function for the Ising 
model with zero field and constant interaction energies can be reformulated as a tractable 
problem for planar graphs (see [TU [28j [22] ). Considerable work has been done investigating 
the computational complexity of the variable field Ising model under under several condi- 
tions, notably by Goldberg and Jerrum J2U], and Jerrum and Sinclair [25]. These include 
complexity classifications and, where possible, approximation algorithms for the Ising model 
with different restrictions on the interaction energies and magnetic fields. In particular there 
is no fully polynomial randomized approximation scheme (FPRAS) in the antiferromagetic 
case, but there is a FPRAS in the ferromagnetic case provided that all the local magnetic 
field values have the same sign. Without this restriction on the magnetic field values, the 
problem again becomes intractable. 



However, Theorems 5.1 and 5.8 imply that the Ising model with constant interaction 
energies and a constant magnetic field vector (without all entries necessarily having the 
same sign), is an evaluation of the [/-polynomial. Thus, by the results of Noble [38], it may 
be computed in polynomial time for graphs with bounded tree width. It is very likely that 
this result may be improved by restricting to the q = 2 case of the Ising model, and also 
likely that these results might be extended to variable interaction energies. Applying the 
theory in the other direction, the complexity results for approximating the Ising model with 
external field from e.g. [201 125] now immediately apply to give computational complexity 
information for the K-polynomial when q = 2. 



6.2. Zeros and phase transitions. In addition to the computational complexity results 
discussed above, there has been considerable research into the zeros of the chromatic and 
multivariate Tutte polynomial. From a graph theory perspective, this was traditionally mo- 
tivated by graph colouring questions, since the chromatic polynomial of a graph G, when 
evaluated at a non-negative integer q, gives the number of ways to properly colour G using q 
colours. From a statistical mechanics perspective, the interest stems from phase transitions 
which may be identified by (accumulation points of) zeros in the partition function. (See, for 
example, J5[ H21 E2 H6[ H7] and the references therein.) For example, considerable effort has 
been devoted to clearing regions of the complex plane, particularly those containing intervals 
of the real axis, of roots of the chromatic polynomial, thus precluding any phase transitions. 
However, these investigations, particularly when theoretical as opposed to computer simula- 
tions, are largely in the absence of an external magnetic field. The authors are not aware of 
any investigations of the zeros of the W-polynomial. It is likely that the techniques devel- 
oped for the classical and multivariable Tutte and chromatic polynomials may be adapted 
to the U-, V- and ly-polynomials and thus extended to Potts models with external fields. 

21 



6.3. Knot theory connections. Connections between knot theory and Potts model were 
first noted over twenty years ago by Jones [26J, and were made concrete by Kauffman [29||30] 
shortly thereafter. Such relations between statistical mechanics and knot theory were soon 
found for other quantum knot invariants by Jones [27j and Turaev [1H], and have since 
been explored by many others. We refer the reader to either Jones' paper [27] or Wu's 
survey article [56] for an overview of connections between statistical mechanical models and 
knot invariants. The relations between statistical mechanical models (particularly the Potts 
model and ice-type models) are of fundamental and continued importance in knot theory. 
For example there has been recent interest in using topological versions of the Potts model 
to explore knot invariants (see [35J and the references therein). 

Each of the relations between physics and knot theory mentioned above uses a graph 
based on the link diagram. The [/-polynomial (and therefore the W- and V-polynomials) 
can be used to construct knot invariants in a very different way. In fact, it is this different 
construction, due to Chmutov, Duzhin and Lando in the sequence of papers [121 H3l E] , that 
led Noble and Welsh to define the U- and jy-polynomials. (The U- and ^-polynomials arose 
through an investigation of the combinatorics behind Chmutov, Duzhin and Lando's work 
on Vassiliev invariants.) This construction of knot invariants via the [/-polynomial uses the 
theory of quantum and Vassiliev knot invariants (see [39], for example, for an introduction 
to quantum and Vassiliev knot theory). In short, by applying the [/-polynomial to the 
intersection graph of a chord diagram one obtains a weight system (a map A — > Q, where 
A is the algebra of chord diagrams from the theory of Vassiliev invariants). One obtains a 
knot invariant by composing this weight system with the Kontsevich invariant. The resulting 
knot invariants were categorized by Lieberum in [32] • By the results presented in this paper, 
these Vassiliev invariants provide a new connection between knot invariants and the Potts 
model. 

It is very provocative to have two very different connections between the Potts model and 
quantum knot invariants. The problem of fully understanding the connections between these 
two applications of statistical mechanics to knot theory is one that the authors consider to 
be important. 

6.4. Other Hamiltonians. While we have given examples of several Potts models that 
may be unified by the V-polynomial, this list is not exhaustive, being only intended to 
illustrate the applications and techniques. Because of the generality of the indexing set, 



even Theorem |4.2| might be adapted to other applications. A more ambitious direction 
would be determining if non-linear terms in the Hamiltonian, such as the squared differences 
that appear in some biological models (see [181 SO] for example), might be assimilated into 
this theory in some way. 

6.5. Combinatorial properties. Although the motivation for defining the V-polynomial 
comes from statistical mechanics, it is an interesting combinatorial object in its own right, 
with potential applications in knot theory as well as statistical mechanics. Thus, further 
study of its combinatorial properties may prove fruitful. The papers [3] and [52] are par- 
ticularly relevant in this regard. [3] discovers a number of graph theoretical properties that 
are encoded by the Ising model with constant interaction energies and magnetic field. Thus, 
all these properties are also encoded by the V-polynomial. A natural direction of investi- 
gation is to extend such results via the V-polynomial to more general situations. The O 
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polynomial of [52] is especially interesting, and since it is equivalent to the RFIM, it is also a 



specialization of the V-polynomial by Theorem 5.8 This relation provides the appropriate 
transform of the polynomial to relate it to the Tutte polynomial, and also gives a deletion- 
contraction reduction for O that holds for non-constant vertex weights. Furthermore, one 
of the results of [52] is that the polynomial gives the ratio of the RFIM to the Bethe 
approximation, an important estimate of the partition function. Thus, since both and 
the RFIM are evaluations of the V-polynomial, the Bethe approximation may be expressed 
as a ratio of V-polynomials. The connection between the V-polynomial and polynomial 
warrants further exploration. 
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